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Abstract 

Defects change essentially not only electronic but also chemical properties of graphene being centers of 

its chemical activity. Their functionalization is a way to modify electronic and crystal structure of graphene 

which may be important for graphene-based nanoelectronics. Using hydrogen as an example, we have 

simulated a chemistry of imperfect graphene for a broad class of defects (Stone- Wales defects, bivacancy, 

nitrogen substitution impurity, and zigzag edges) by density functional calculations. We have studied also 

^^' ■ an effect of finite width of graphene nanoribbons on their chemical properties. It is shown that magnetism 

f^ ' at graphene edges is fragile with respect to oxidation and, therefore, a chemical protection of the graphene 

rvQ , edges may be required for application of graphene in spintronics. At the same time, hydrogenation of the 

.| ' Stone- Wales defects may be a perspective way to create magnetic carbon. 

^ ' Graphene [H 13 H] is a novel and very perspective material for nanoelectronics [H H] . Effect of various 

imperfections of the crystal lattice, such as Stone- Wales defects (SW) [5], vacancies [5], substitution atoms [71IH|, 
etc., on physical properties of graphene was studied theoretically [101 HH [El [E] . Recently, first experimental 
observations of defects in graphene have been reported [HI [151 [16] . Earlier, the defects in carbon nanotubes have 
been intensively investigated [T71 [TSl [IHl HOI HI] ■ It is known that defects can influence essentially on chemistry 
C/3 , of graphene, and there are several theoretical considerations of chemical functionalization of imperfect graphene 

[12 im US] ; also, the defects were experimentally used as attractors for impurities [E] . 

It is still unclear, however, what are typical defects in real graphene samples and what is their concentration. 
Further studies would be important both scientifically and technologically, to find a way to produce defect-free 
graphene or graphene with desired defects to modify electronic properties, like in conventional semiconductors. 
C^ ' The chemical functionalization is probably the most suitable way to detect imperfections in graphene crystal 

lattice. However, only chemisorption of individual atoms or their pairs were considered before, whereas it was 
demonstrated that, say, for the case of hydrogenation of ideal graphene a complete coverage turned out to be 
the energetically most favorable so, the most probably, impurities would tend to form the whole clusters on the 
Q ' top of imperfect graphene [231 [5S] . Therefore a systematic study of chemisorption on defects in graphene for a 

CJ ! broad range of coverage is desirable; this is the subject of the present work. 

Edges of graphene sheets which are broadly discussed now, due to their potential relevance for field- 

^vg I effect transistors and spintronic devices [2S1 [HI HHl [2S [301 [SH 132 can be also considered as defects. Their 

^ ' functionalization and its effect on the electronic structure were discussed in a series of theoretical works 

in! [3Tl[33l|34l|35l[36l[33[38l[39l[3ni|41], chemical activity of the graphene edges being discussed [42], but only 

in ■ the chemical functionalization of the edges themselves were considered rather than of the whole graphene sheet 

with the edges. Actually, the edges change chemical activity of graphene in an extended region, as will be 

demonstrated here. This is important to study stability of magnetism at the edges and to evaluate a width of 

a "near-edge layer" where properties of graphene are essentially different from those in the bulk. 

We will use hydrogenation as an example of the chemical functionalization of defects in graphene. For the 
f ^ I case of ideal graphene general factors determining the chemisorption were first studied for the case of hydrogen 

^^ and appeared to be applicable also in other situations, such as graphene oxide [53]. One can expect therefore 
that, at least, some of peculiarities which are found here can be of a more general relevance. 
^^ ■ We used the pseudopotential density functional SIESTA package for electronic structure calculations [331 [33] 

H I with the generalized gradient approximation for the density functional [45], with energy mesh cutoff 400 Ry, 

and /c-point 8x8x2 mesh in Monkhorst-Park scheme [36] for modeling of defects, and 2x12x2 for modeling of 
nanoribbons. During the optimization, the electronic ground state was found self-consistently by using norm- 
conserving pseudopotentials for cores and a double-(^ plus polarization basis of localized orbitals for carbon and 
oxygen, and double-^ one for hydrogen. Optimization of the bond lengths and total energies was performed 
with an accuracy 0.04 eV /A and 1 meV, respectively. This method is frequently used for computations of 
electronic structure of graphene [9] [28l [29l [231 [25] . When drawing the pictures of density of states, a smearing 
by 0.2 eV was used. 

Chemisorption energies were calculated by standard formulas used, e.g., earlier for the cases of chemisorption 
of hydrogen j24j oxygen and hydroxyl groups ,25j on graphene and solution of carbon in 7-iron ^47j . Thus, energy 
of the chemisorption of single hydrogen atom at graphene with Stone- Wales defect (Fig. [T^) was calculated as 



E/orm=EsvK+ff-Esw-E_f/2/2 where Esw+h was the total energy of the supercell with chemisorbed hydrogen 
found by self-consistent calculations after optimization of geometric structure, Fjsw was the total energy of 
the graphene supercell with the Stone- Wales defect found in a similar way, and Eh2 is the energy of hydrogen 
molecule. The latter choice describes the most adequately a process of spontaneous hydrogen dissociation and, 
thus, stability of graphene with respect to the hydrogenation. That is why this definition is broadly used when 
discussing hydrogen storage [IHIHH]- In principle, a definition of the hydrogenation energy depends on what 
specific compound is used for the hydrogenation. This may shift the chcmisorption energy values of hydrogen in 
Figs, m HI m and 1101 For example, for the hydrogenation by hydrogen plasma all formation energies presented 
below will be lower by 104.206 kcal/mol, or 2.255 eV per atom, which is the experimental bond dissociation 
energy for H2. 

We used graphene supercell containing 50 carbon atoms, similar to that used in our previous work [24| . 
It was demonstrated there that the chcmisorption energy for two hydrogen atoms is close to that for two 
independent single atoms if the distance between the atoms is larger than 1 nm (approximately, four graphene 
lattice constants). Similar results have been obtained for the case of hydroxy groups ^5\. This means that a 
supercell provided the distance between defects larger than 1 nm is sufficient to model the chcmisorption on 
individual defects. 

Optimized geometric structures of graphene supercells with various kinds of defects are presented in Fig. 
[TJ To simulate graphene nanoribbons we used the stripes with 22, 44, and 66 carbon atoms per cell, with the 
width 2.20, 4.54, and 6.88 nm, respectively. 

We start with the calculations of activation energy, that is, the chcmisorption of a single hydrogen atom. 
For pure graphene this energy is 1.5 eV [21] and for all kind of defects under consideration it turns out to be 
lower, namely, 0.30 eV for SW, 0.93 eV for bivacancy, and 0.36 eV for substitution impurity of nitrogen (in the 
latter case we mean the energy of chcmisorption on a carbon neighbor of the nitrogen impurity) . 

We have calculated also the chcmisorption energy for different nonequivalent atoms of the lattice containing 
the SW defect (Fig. [^h-)- One can see that the chcmisorption energy at the SW defect is minimal. At the 
same time, in a whole area surrounding the defect the chcmisorption energy is lower than in perfect graphene 
[53] . For distant enough atoms the energy is very close to the values found in Ref. [23] which confirms that the 
supercell with 50 atoms is sufficient to describe the chcmisorption on the single defect. 

The computational results mean that defects in graphene are centers of chemical activity. As a next step, 
we investigate an interaction of a pair of hydrogen atoms with defects (Fig. |3^). Again, the defects decrease the 
chcmisorption energy of the pair in comparison with the case of ideal graphene. We have calculated also the pair 
formation energy for the case when the first hydrogen atom is chemisorbed on the SW defect and the second 
one - on one of the surrounding carbon atoms (Fig. |2)d) . The chcmisorption energies are presented there for the 
case of one-side chcmisorption; for the first and second neighbors of the SW defect the one-side chcmisorption is 
energetically more favorable [T21I2HI23] whereas for farther atoms the difference with the two-side chcmisorption 
energies becomes negligible [23] • Similar to the case of single hydrogen atom discussed above a region is formed 
around the defect with a diameter of order of 1 nm where the chcmisorption energies are lower than in the bulk. 

These results may be essential for the problem of magnetism in graphene. Chcmisorption of a single hydrogen 
atom on graphene results in appearance of magnetic moments [12l |50l |5T] whereas close hydrogen pairs on 
graphene turn out to be nonmagnetic, the pairs being more energetically favorable than independent atoms 
[241 1521 I53j . At the same time, there is an energy barrier which two distant hydrogen atoms on graphene should 
overcome to form the pair [S3j . If one of these atoms are bonded with the SW defect the barrier becomes higher 
by 0.44 eV which makes the pair formation at room temperature rather difficult. Thus, the hydrogenation of 
the SW defects may be a perspective way to create magnetic graphene. It can be done, for example, using 
reactions with superacids, similar to the case of fullerenes |55) . 

For SW and bivacancies the chcmisorption energy of single hydrogen atom is negative which makes these 
defects centers of stable chemisorption of hydrogen. At further increase of the number of hydrogen atoms the 
chcmisorption energy decreases reaching a minimum for the most stable configuration (for the case of SW it is 
shown in Fig. [Sh), then increases up to a local maximum (the corresponding structure is shown in Fig. [3j:) and 
then decreases again, with the global minimum for complete coverage. Contrary, in the case of ideal graphene 
the chemisorption energy decreases monotonously with the hydrogen concentration and clusters of chemisorbed 
hydrogen are not stable [b& . Interestingly, for the case of complete coverage the binding energy is smaller for 
the hydrogen on graphene with defects than on the perfect one. This means that completely hydrogenated 
graphene (graphane [231ISI]) is less stable with defects than without them. 

Our results show that, whereas at the hydrogenation of pure graphene there is single potential barrier 
to overcome, that is, the adsorption of the first atom, for imperfect graphene several barriers can exist (see 
Fig. 131). The first one, again, corresponds to the activation energy and the other ones - to the transition 



from hydrogenation of the region around defect with the minimal chemisorption energy (see Fig. [2]) to the 
hydrogenation of the whole graphene sheet. Studying experimentally kinetics of this process one can extract an 
information about presence of defects in graphene and their type. At the same time, graphene with defects can 
be a good model of activated carbon [511 [SH] ■ The fact that the adsorption and desorption energies are smaller 
than 0.5 eV per hydrogen atom makes graphene with defects a perspective catalyst. 

Let us consider now the hydrogenation of graphene edges. We focus on the zigzag edges which are especially 
interesting due to their possible half-metallic ferromagnetic state [281 ESI |30] . For the nanoribbons of width 
2.20, 4.54, and 6.88 nm, the chemisorption energies for a single hydrogen atom (Fig. O)) are found to be -2.39, 
-2.76, and -2.82 eV, respectively. High chemical activity of the edges, in comparison with the bulk graphene [25 
is explained by a presence of unpaired electrons due to broken bonds (Fig. [5^) . Only one unpaired electron is 
involved in the bond with the single hydrogen atom whereas the second one can participates in the second bond 
(Fig. [5t). The chemisorption of the second hydrogen atom leads to strong local distortions of graphene, due 
to a change of hybridization from sp^ (favorable for a planar geometry) to sp'^ (favorable for three-dimensional 
structure like in diamond). These distortions require some energy which is higher than the energy gain due 
to pairing of all electrons and, as a result, an energy barrier about 1.5 cV is formed for all nanoribbons under 
investigation. The cohesive energy of pure graphene is rather strongly dependent on the ribbon width (Fig. 
[6]). This result in a strong dependence of the chemisorption energy on the ribbon width as well. The energy of 
chemisorption of two hydrogen atoms at the edges is negative, in contrast with the case of perfect bulk graphene 
[53] which makes the edges another natural centers of chemical activity. Similarly, zigzag edges can adsorb two 
fluorines [61] . or two hydroxyl groups, or a pair of hydroxyl group and hydrogen, the latter case being related 
to the decomposition of a water molecule (Fig. [7]). The formation energy per H2O molecule corresponding to 
this process is -2.42 eV (241 kJ/mol) for the ribbon width 2.20 nm which makes the decomposition of the water 
more favorable than the chemisorption of single hydrogen atom. 

In a case of chemisorption of oxygen atoms, both unpaired electrons are involved but the hybridization 
remains sp^ and local geometry therefore remains flat (Fig. [SJi) . The chemisorption energy for the oxygen atom 
is much lower than for hydrogen or hydroxyl; for the ribbon width 2.2 nm it turns out to be -4 eV. Recently, 
peaks corresponding to C=0 chemical bonds have been found in the core-level spectra of carbon in graphene 
nanoribbons |60j . They were related to edge groups, based on ratio of intensities of these peaks and of the main 
peaks corresponding to sp^ hybridization. 

Since all unpaired electrons are bonded at the complete oxygenation of the edge, the magnetism disappears 
and electronic structure becomes purely metallic (see Fig. |S]) , instead of half-metallic ferromagnetic for undoped 
zigzag edges ^8\. Note that the decomposition of the water molecule described above will also suppress the 
magnetism of the zigzag edges, as was observed at fluorination of the graphene edges [6T] . 

Now we consider further hydrogenation of graphene, from the edges to bulk. The computational results 
presented in Fig. [5] show that the chemisorption of hydrogen is the most energetically favorable near the edges 
and less favorable in the middle of the nanoribbon where the formation energy is close to that for the bulk 
(1.44 eV [Mj)- The chemisorption energy is very sensitive to the distortions created by adsorbed hydrogen 
atoms. It was shown earlier for the bulk graphene |24j that the difference in chemisorption energies for different 
configurations of a pair of hydrogen atoms can be as large as 0.45 eV. In graphene nanoribbons these effects 
seem to be even stronger. In particular, geometric frustrations due to the atomic distortions lead to a zigzag 
dependence of the chemisorption energy on atomic position (Fig. [5]). At the chemisorption of single hydrogen 
atom interactions between the hydrogen- induced distortions and those due to edge itself is important. This 
effect decays when moving from the edge at distances of order of 1 nm. As well as in the case of SW defect, the 
pair formation can be more difficult for the nanoribbon than for bulk graphene, due to larger barrier energy. 
This may be essential, again, for the problem of hydrogen-induced magnetism in graphene. 

The results are shown in Fig. (TUl One can see that, irrespective to the ribbon width, at first the formation 
energy grows by steps which corresponds to a preferable hydrogenation of graphene by pairs of atoms 23). 
This growth continues at the depth of 8 carbon atoms (0.85 nm, which is slightly larger than for the single 
atom due to two-side hydrogenation) and then the energy increases approximately linearly with the coverage, 
as in the case of infinite graphene sheet [24] . This means that the hydrogenation of graphene proceeds in three 
steps, namely, the hydrogenation of the edge, then, of the near-edge layer with a width of order 1 nm, and, 
at last, of the bulk. The chemical functionalization of the near-edge layer effects essentially on the electronic 
structure of graphene nanoribbons (see Fig. [S]) . It is worth to note that for the ribbons under consideration the 
hydrogenation energy is negative which makes them chemically active in their whole area. This can infiuence 
the electronic structure as was mentioned earlier [26] . 

To conclude, let us summarize the main results. Defects in graphene crystal lattice, as well as the edges 
of graphene, change drastically a scenario of its chemical functionalization. For the ideal graphene, only sin- 



gle energy barrier related with the chemisorption of the first atom is relevant and complete functionalization 
corresponds to the global energy minimum. At the same time, for "realistic" graphene with imperfections and 
edges local energy minima are formed corresponding to local functionalization of a region with radius of order 
half of nanometer and further functionalization is stopped by presence of another energy barrier. The local 
energy minima make migration of chemisorbed atoms and groups more difficult. In the cases of hydrogen and 
other atoms and groups with single unpaired electron, this can stabilize a magnetic state preventing formation 
of nonmagnetic pairs. We have found that magnetic state of edges in graphene nanoribbons may be unstable 
with respect to oxidation and water dissociation at the edges. Therefore, keeping in mind potential applications 
of this magnetism in spintronics special methods of chemical protection of the edges should be developed. 
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Figure 1: Optimized geometric structures for graphene supercell containing the Stone-Wales defect (a), bivacancy (b), 
and nitrogen substitution impurity (c). 
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Figure 2: Formation energies (in eV per hydrogen atom) for nonequivalent sites (red circles) of graphene lattice with 
the SW defect for chemisorption of the first (a) and second (b) hydrogen atoms. 
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Figure 3: Optimized geometric structures for the Stone-Wales defect functionalized by 2 (a), 6 (b), 14 (c) hydrogen 
atoms and completely covered by hydrogen (d). Carbon and hydrogen atoms are shown by green and violet circles, 
respectively. 
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Figure 4: Energy formation as a function of coverage for grapehne sheet containing the Stone- Wales defect (a), bivacancy 
(b), and nitrogen substitution impurity (c). The blue dashed line represents the results for the ideal infinite graphene 
sheet. 
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Figure 5: A sketch of functionalization of the zigzag edges of graphene; (a) Initial state of edges, unpaired electrons 
are shown in blue; (b), (c), and (d) the edge functionalized by a single carbon atom (yellow), by pair of hydrogen atoms, 
and by oxygen atom (red), respectively. 
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Figure 6: Total energy difference between infinite graphene sheet and graphene nanoribbon as a function of the 
nanoribbon width. 




Figure 7: Optimized geometry of the graphene nanoribbon of width 2.20 nm with the edges functionalized by hydrogen 
atoms and hydroxy groups resulted from water decomposition. 
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Figure 8: Densities of states for graphene nanoribbon of width of 2.2 nm with its edges functionalized by: (a) oxygen 
atoms and (b) hydrogen atoms; (c) the case of complete coverage of six first layers from the edges On insets optimized 
geometric structures of fuctionalized nanoribbons corresponding with dencities of states. 
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Figure 9: Chemisorption energy of hydrogen as a function of distance d from the edge. Inset: optimized geometric 
structure of the system under consideration. 
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Figure 10: Formation energy as a function of coverage of graphene nanoribbons by hydrogen as a function of the 
nanoribbon width (red sohd hue). The blue dashed line presents the results for an ideal infinite graphene sheet. 
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